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q 1 The strong CP problem 

<d : 

The axion was postulated nearly two decades ago |]J to explain why the strong interactions 
conserve P and CP in spite of the fact that the weak interactions violate those symmetries. 
^ ■ Consider the Lagrangian of QCD: 

& • 

1 n Bo 2 

Cqcd = — A Gl v G a ^ + £ [Vj-fiDM - {m jq %q Rj + h.c.)] + ^G%G^ . (1.1) 

i=i 

The last term is a 4-divergence and hence does not contribute in perturbation theory. That 
term does however contribute through non-perturbative effects associated with QCD 
instantons 0. Such effects can make the physics of QCD depend upon the value of 9. 
Using the Adler-Bell- Jackiw anomaly [|J , one can show that 9 dependence must be there if 
none of the current quark masses vanishes. (If this 9 dependence were absent, QCD would 
have a Ua{X) symmetry and would predict the mass of the r/' pseudo-scalar meson to be 
less than ^JZm^ w 240 MeV ||, contrary to observation.) One can further show that QCD 
depends upon 9 only through the combination of parameters: 

9 = 9 - arg(rai, m 2 , . . .m n ) (1.2) 
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If 9 ^ 0, QCD violates P and CP. The absence of P and CP violations in the strong 
interactions therefore places an upper limit upon 9. The best constraint follows from the 
experimental bound |J on the neutron electric dipole moment which yields: 9 < 10~ 9 . 

The question then arises: why is 9 so small? In the standard model of particle inter- 
actions, the quark masses originate in the electroweak sector of the theory. This sector 
must violate P and CP to produce the correct weak interaction phenomenology. There is 
no reason in the standard model to expect the overall phase of the quark mass matrix to 
exactly match the value of 9 from the QCD sector so that 9 < 10 -9 . In particular, if CP 
violation is introduced in the manner of Kobayashi and Maskawa , the Yukawa couplings 
that give masses to the quarks are arbitrary complex numbers and hence arg det m q and 
9 have no reason to take on any special value at all. 

The problem why 9 < 10~ 9 is usually referred to as the "strong CP problem". The 
existence of an axion would solve this problem. But, before we talk about axion physics, 
let's mention that there are other solutions. Setting m u = removes the ^-dependence of 
QCD and hence the strong CP problem as well. However, m u = may cause problems 
with the successful current algebra relations among pseudo-scalar meson masses. I refer the 
reader to refs. || |9| for recent discussions of the issues involved. Another type of solution 
involves the assumption that CP and/or P is spontaneously broken but is otherwise a 
good symmetry. In this case, 9 is calculable and may be arranged to be small [ 10| . Finally, 



let's emphasize that the strong CP problem need not be solved in the low energy theory. 
Indeed, as Ellis and Gaillard |Tl|] pointed out, if in the standard model 9 = near the 



Planck scale then 9 <C 10 9 at the QCD scale. 

Peccei and Quinn JT2J proposed to solve the strong CP problem by postulating the 



existence of a global Upq(1) quasi-symmetry. To do its job, Upq(1) must be a symmetry 
of the theory at the classical (i.e., at the Lagrangian) level, it must be broken explicitly by 
those non-perturbative effects that make the physics of QCD depend upon 9, and finally 
it must be spontaneously broken. The axion [13| is the pseudo-Nambu-Goldstone boson 
associated with the spontaneous breakdown of Upq(1). One can show that, if a Upq(1) 
quasi-symmetry is present, then 

at x) 

9 = 9 - arg{m x . . .m n ) — , (1.3) 

la 

where a(x) is the axion field and f a , called the axion decay constant, is of order the vacuum 
expectation value (VEV) which spontaneously breaks Upq(1). It can further be shown 
[TJJ] that the non-perturbative effects that make QCD depend upon 9 produce an effective 



potential V(9) whose minimum is at 9 = 0. Thus, by postulating an axion, 9 is allowed 
to relax to zero dynamically and the strong CP problem is solved. 



The properties of the axion can be derived using the methods of current algebra (15 
The axion mass is given in terms of f a by 

m a ~ 0.6 eV W7( ^ eV . (1.4) 

J a 

All the axion couplings are inversely proportional to f a . For example, the axion coupling 
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to two photons is: 

C arr = ■ B (1.5) 
where and B are the electric and magnetic fields, a is the fine structure constant, and g y 



is a mo del- dependent coefficient of order one. g y = 0.36 in the DFSZ model |L6| whereas 



<7 7 = —0.97 in the KSVZ model |T7[ . A priori the value of / a , and hence that of m a , is 
arbitrary. However, searches for the axion in high energy and nuclear physics experiments 
combined with astrophysical constraints, the latter derived by considering the effect of 
the axion upon the lifetimes of red giants and SN1987a, rule out m a > 10~ 3 eV [0. In 
addition, as will be discussed in section II, cosmology places a lower limit on m a of order 
10~ 6 eV by requiring that axions do not overdose the universe. 



2 Dark matter axions 

For small masses, axion production in the early universe is dominated by a novel mechanism 
|T8|j . The crucial point is that the non-perturbative QCD effects that produce the effective 



potential V{9) are strongly suppressed at temperatures high compared to Aqcd 101 • At 
these high temperatures, the axion is massless and all values of (a(x)) are equally likely. 
At T ~ 1 GeV, the potential V turns on and the axion field starts to oscillate about a 
CP conserving minimum of V. These oscillations do not dissipate into other forms of 
energy because, in the relevant mass range, the axion is too weakly coupled for that to 
happen. The oscillations of the axion field may be described as a fluid of axions. The 
typical momentum of the axions in the fluid is the inverse of the correlation length of the 
axion field. Because that correlation length is of order the horizon, we have p a ~ (10~ 6 
sec)" 1 ~ 10 -9 eV at T ~ 1 GeV, and p a ~ -R" 1 afterwards. R is the cosmological scale 
factor here. Thus the axion fluid is very cold compared to the ambient temperature. 

Let me briefly indicate how the present cosmological energy density of this axion fluid 
is estimated. Let ^p(x) be the complex scalar field whose VEV v spontaneously breaks 
Upq(1). At extremely high temperatures, the Upq(1) symmetry is restored. It becomes 
spontaneously broken when the temperature drops below a critical value Tpq of order v. 
Below Tpq, the axion field a(x) appears as the phase of the VEV of </?: (ip(x)) = ve m ^' v . 
We must now distinguish two cases. Either inflation occurs with reheat temperature 
below T P Q, or not (i.e., inflation does not occur or it occurs with reheat temperature 
above Tpq). In the first case, inflation homogenizes the axion field and there is only 
one contribution to the axion cosmological energy density, the contribution from so-called 
"vacuum misalignment" . In the second case, there are additional contributions from axion 
string and axion domain wall decay. Only the contribution from vacuum misalignment is 
discussed in any detail here. 

When the axion mass turns on near the QCD phase transition, the axion field starts to 
oscillate about one of the CP conserving minima of the effective potential. The oscillation 
begins approximately at cosmological time t± such that tim a (T(ti)) = 0(1), where m a (T) is 
the temperature dependent axion mass. Soon after t\, the axion mass changes sufficiently 
slowly that the total number of axions in the oscillations of the axion field is an adiabatic 
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invariant. T\ = T(t\) has been estimated to be of order 1 GeV. The number density of 
axions at time ti is 

n^O^m^^i))^! (2.1) 

where / = jr is the axion decay constant introduced earlier. iV is an integer which 
expresses the color anomaly of Upq(1). N also equals the number of CP conserving 
vacua [EjJ at the bottom of the 'Mexican hat' potential, i.e., in the interval < - < 2tt. 
In Eq. (2.1), we used the fact that the axion field a(x) is approximately homogeneous on 
the horizon scale t\. Wiggles in a(x) which entered the horizon long before t± have been 
red-shifted away [21]. We also used the fact that the initial departure of a(x) from the 
nearest minimum is of order = f a . The axions of Eq. (2.1) are decoupled and non- 
relativistic. Assuming that the ratio of the axion number density to the entropy density 



is constant from time t\ till today, one finds [18 



/0.6 HT 5 eV\^ /200 MeV\' /75 km/s ■ Mpc\ 2 
a ~ [ m a J [ Aqcd J [ H~o ) ( ' } 

for the ratio of the axion energy density to the critical density for closing the universe. H 
is the present Hubble rate. Eq. (2.2) implies the bound m a > 10~ 6 eV. 

However, it should be emphasized that there are many sources of uncertainty in the 
estimate of Eq. (2.2). The axion energy density may be diluted by the entropy release 



from heavy particles which decouple before the QCD epoch but decay afterwards ||22|| , 
or by the entropy release associated with a first order QCD phase transition. On the 
other hand, if the QCD phase transition is first order [p^| , an abrupt change of the axion 
mass at the transition may increase Q a . If inflation occurs with reheat temperature less 
than Tpq, there may be an accidental suppression of Q a because the homogenized axion 
field happens to lie close to a CP conserving minimum. Because the RHS of Eq. (2.2) is 
multiplied in this case by a factor of order the square of the initial vacuum misalignment 
angle ^p-N which is randomly chosen between — it and +7r, the probability that Vl a is 
suppressed by a factor x is of order y/x. This rule cannot be extended to arbitrarily small 
x however because quantum mechanical fluctuations in the axion field during the epoch of 
inflation do not allow the suppression to be perfect |p4|| . If inflation occurs with reheating 



temperature larger than Tpq or if there is no inflation, there are contributions to Q a from 



axion string [£5] and axion domain wall decay in addition to the contribution, Eq. (2.2), 
from vacuum misalignment. The author and his collaborators [EEJ have estimated each of 
these additional contributions to be of the same order of magnitude as that from vacuum 



misalignment. Others [BR |27j have estimated that the contribution from axion string 



decay dominates over that from vacuum misalignment by a factor 100 or a factor 10. 

The axions produced when the axion mass turns on during the QCD phase transition 
are cold dark matter (CDM) because the axions are non-relativistic from the moment of 
their first appearance at 1 GeV temperature. Studies of large scale structure formation 



support the view that the dominant fraction of dark matter is CDM ||28|| . Moreover, any 
form of CDM necessarily contributes to galactic halos by falling into the gravitational wells 
of galaxies. Hence, there is excellent motivation to look for CDM candidates as constituent 
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particles of our galactic halo, even after some fraction of our halo has been demonstrated 



to be in MACHOs 29 or some other form. 



Finally, let's mention that there is a particular kind of dumpiness [30] which affects 
axion dark matter if there is no inflation after the Peccei-Quinn phase transition. This 
is due to the fact that the dark matter axions are inhomogeneous with 5p/p ~ 1 over 
the horizon scale at temperature T\ ~ 1 GeV, when they are produced at the start of 
the QCD phase-transition, combined with the fact that their velocities are so small that 
they do not erase these inhomogeneities by free-streaming before the time t eq of equality 
between the matter and radiation energy densities when matter perturbations can start 
to grow. These particular inhomogeneities in the axion dark matter are immediately in 
the non-linear regime after time t eq and thus form clumps, called 'axion mini-clusters' |30[ . 
mc ~10- 13 M Q 



These have mass M mc ~ 10 13 M M and size l mr ~ 10 12 cm. 



3 The cavity detector of galactic halo axions 

Axions can be detected by stimulating their conversion to photons in a strong magnetic field 



31]. The relevant coupling is given in Eq. (1.5). In particular, an electromagnetic cavity 
permeated by a strong static magnetic field can be used to detect galactic halo axions. The 
latter have velocities (3 of order 10~ 3 and hence their energies E a = m a + |m a /5 2 have a 
spread of order 10~ 6 above the axion mass. When the frequency u = 2k f of a cavity mode 
equals m a , galactic halo axions convert resonantly into quanta of excitation (photons) of 
that cavity mode. The power from axion — > photon conversion on resonance is found to 



be 31 



VB 2 oPa C—Mm(Q L ,Q a ) 
* U) rn a 



0.5 10" 26 Watt 



V \ ( B ° \c( 9 ^ 
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■ (ri5%)(2«) Mi ^«»' (") 

where V is the volume of the cavity, B is the magnetic field strength, Ql is its loaded 
quality factor, Q a = 10 6 is the 'quality factor' of the galactic halo axion signal (i.e. the 
ratio of their energy to their energy spread), p a is the density of galactic halo axions on 
Earth, and C is a mode dependent form factor given by 



C 



. 2 



BlV f v d 3 xe\E w \ 2 



(3.2) 



where B (x) is the static magnetic field, Eu{x)e l03t is the oscillating electric field and e is 
the dielectric constant. 

Because the axion mass is only known in order of magnitude at best, the cavity must 
be tunable and a large range of frequencies must be explored seeking a signal. The cavity 
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Figure 1: Regions in mass - (coupling/mass) 2 space which have been ruled out by the 
RBF and UF experiments (hatched) and which the LLNL experiment expects being able 
to rule out (shaded border). The relation between coupling and mass in the DFSZ and 
KSVZ models is also shown. 

can be tuned by moving a dielectric rod or metal post inside it. Using Eq. (3.1), one finds 
that to perform a search with signal to noise ratio s/n, the scanning rate is: 



df 12GHz (An\ 
dt 




(3.3) 

where T n is the sum of the physical temperature of the cavity plus the noise temperature of 
the microwave receiver that detects the photons from a — > 7 conversion. Eq. (3.3) assumes 
that Ql < Qa and that some strategies have been followed which optimize the search rate. 
The best quality factors attainable at present, using oxygen free copper, are of order 10 5 
in the GHz range. To make the cavity of superconducting material is probably not useful 
since it is permeated by a strong magnetic field in the experiment. 

Eq. (3.3) shows that a galactic halo search at the required sensitivity is feasible with 
presently available technology, provided the form factor C can be kept at values of order one 
for a wide range of frequencies. For a cylindrical cavity and a longitudinal magnetic field, 
C = 0.69 for the lowest TM mode. The form factors of the other modes are much smaller. 
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The resonant frequency of the lowest TM mode of a cylindrical cavity is /=115 MHz (jjf) 
where R is the radius of the cavity Since 10" 6 eV = 2tt (242 MHz), a large cylindrical 
cavity is convenient for searching the low frequency end of the range of interest. To extend 
the search to high frequencies without sacrifice in volume, one may power-combine many 
identical cavities which fill up the available volume inside a magnet's bore P3"| , ! SA\. This 
method allows one to maintain C = 0(1) at high frequencies, albeit at the cost of increasing 
engineering complexity as the frequency, and hence the number of cavities, is increased. 

Pilot experiments were carried out at Brookhaven National Laboratory |35|] and at the 
University of Florida [pq| . The (magnetic field) 2 x volume provided by the magnets used 
in these experiments were relatively low: BqV = 0.36T 2 m 3 and 0.45T 2 m 3 respectively for 
Brookhaven (RBF) and Florida (UF). Fig. I shows the limits that these experiments placed 
on the square of the coupling g ail = ^j- as a function of the axion mass m a assuming 
that the local density of galactic halo axions is p a — \ ■ 10~ 24 gr/cm 3 . The figure also shows 
function of m a in the KSVZ and DFSZ models. 

Second generation experiments are presently under way at Lawrence Livermore Na- 
tional Laboratory (LLNL) |J7] and at Kyoto University [3S|. The LLNL experiment uses 



a much larger magnet {BqV = 12T 2 m 3 ) than the pilot experiments. It also improves the 
noise temperature (T n = 3K vs. T n = 5K for the Florida experiment) although it uses the 
same microwave receiver technology (HEMT amplifiers). The LLNL experiment is also 
the first to use multiple cavity arrays to expand widely the mass range searched. It will 
cover 1.3 < m a < 13/ieV at a level of sensitivity sufficient to discover KSVZ axions if they 
are the constituents of our galactic halo; see Fig. I. The LLNL experiment started taking 
data in Feb. '96 and will run for about three years to cover the stated range. 

The Kyoto experiment has a magnet of size similar to that of the pilot experiments but 
uses a beam of Rydberg atoms to count the photons from a — » 7 conversion. The a —>■ 7 
conversion part is the same as in the other experiments. Single photon counting constitutes 
a dramatic improvement in microwave detection sensitivity. With HEMT amplifiers one 
needs to have thousands of a — > 7 conversions per second and integrate for about 100 sec 
to find a signal in the noise. With single photon counting, a few a — > 7 conversions suffice 
in principle. To build a beam of Rydberg atoms capable of single photon counting is a 
considerable achievement. In addition, the cavity will be cooled by a dilution refrigerator 
down to a temperature (~ 10 mK) where the thermal photon background is negligible. 
The Kyoto experiment will first search near m a = 10~ 5 eV. Its projected sensitivity is 
sufficient to discover DFSZ axions even if their local density is only | of the local halo 
density. 



4 The phase space structure of cold dark matter halos 

If a signal is found in the cavity detector of galactic halo axions, it will be possible to 
measure their energy spectrum with great precision and resolution because all the time 
previously used in searching for the signal can now be used to accumulate data. Hence 
there is good motivation to ask what can be learned about our galaxy from analyzing such 
a signal. 
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E n = (v n /300 km s" 1 )' 



Figure 2: The spectrum of velocity peaks in a typical case (e 
studied in ref. EOf. f n and E n are defined in the text. 



0.2, h = 0.7 and j = 0.2) 



In many past discussions of dark matter detection on Earth, it has been assumed that 
the dark matter particles have an isothermal distribution. Thermalization has been argued 
to be the result of a period of "violent relaxation" following the collapse of the protogalaxy. 
If it is strictly true that the velocity distribution of dark matter particles is isothermal, 
which seems to be a strong assumption, then the only information that can be gained 
from its observation is the corresponding virial velocity and our own velocity relative to 
its standard of rest. If, on the other hand, thermalization is incomplete, a signal in a dark 
matter detector may yield additional information. 

J.R. lpser and I discussed |39|] the extent to which the phase-space distribution of 
cold dark matter particles is thermalized in a galactic halo and concluded that there are 
substantial deviations from a thermal distribution in that the highest energy particles have 
discrete values of velocity. There is one velocity peak on Earth due to dark matter particles 
falling onto the galaxy for the first time, one peak due to particles falling out of the galaxy 
for the first time, one peak due to particles falling into the galaxy for the second time, etc. 
The peaks due to particles that have fallen in and out of the galaxy a large number of 
times in the past are washed out because of scattering in the gravitational wells of stars, 
globular clusters and large molecular clouds. But the peaks due to particles which have 
fallen in and out of the galaxy only a small number of times in the past are not washed 
out. 

I. Tkachev, Y. Wang and I have used the secondary infall model of galactic halo for- 
mation to estimate the local densities and the velocity magnitudes of the dark matter 
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particles in the velocity peaks [^(J. We generalized the existing version of that model to 
take account of the angular momentum of the dark matter particles. In the absence of 
angular momentum, the model produces flat rotation curves for a large range of values 
of a parameter e whose average value may be inferred from the spectrum of primordial 
density perturbations. We find that the presence of angular momentum produces an ef- 
fective core radius, i.e., it makes the contribution of the halo to the rotation curve go to 
zero at zero radius. The model provides a detailed description of the large scale properties 
of galactic halos including their density profiles, their extent and their total mass. Fig. II 
shows the predictions of the model for the average density fractions /„ = p n j p and the 
kinetic energies per unit mass E n of the particles in the highest energy peaks for the case 
where e = 0.2, H = 70 km/sec-Mpc and the average amount of angular momentum, in the 
dimensionless units defined in ref. ||4Q|| , is j = 0.2. The density fractions f n are averages 
over all locations at the same distance (8.5 kpc) from the galactic center as we are. The 
E n are measured in a frame of reference which is not co-rotating with the disk. 
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